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1 Introduction 


The motivation for this paper is to consider the boundary value problem 

— Ay = f in n, y = u on T := dO., (1-1) 

with right hand side / G H~^{yL) and boundary data u G Lp'iT). We assume C to 
be a bounded polygonal domain with boundary T. Such problems arise in optimal control 
when the Dirichlet boundary control is considered in Lp‘iT) only, see for example the 
papers by Deckelnick, Gunther, and Hinze, [11], French and King, [12], May, Rannacher, 
and Vexler, [17], and Apel, Mateos, Pfefferer, and Rosch, [1], On the continuous level 
we even admit more irregular data. 

In Section 2 we analyze several ways how to understand the solution of the boundary 
value problem (1.1) for which we cannot expect a weak solution y G 77^(0). The most 
popular method to solve problem (1.1) is the transposition method that goes back to 
Lions and Magenes [16] and that is based on the use of some integration by parts. This 
formally leads to the very weak formulation: Find y €Y such that 

{y,Av)n = {u,dnv)r-if,v)n Vu G K (1.2) 

with {w,v)g '■= SqWV denoting the L‘^{G) scalar product or an appropriate duality 
product. The main issue is to find the appropriate trial space Y and test space V. In 
the convex case, it turns out that a good choice is K = L^(n) and V = H^(Q) n 77g (Q). 
This case is investigated by many authors, some of them are mentioned in Subsection 
2.1, but we will see that such a choice is not appropriate in the non-convex case (loss 
of uniqueness) and present some remedies (enlarged test spaces) in Subsections 2.2 and 
2.3. We will further show in Subsection 2.2 that in the correct setting the very weak 
solution corresponds to the one obtained by an integral equation technique and hence 
has 77^/^ (n) regularity. 

The main drawback of the very weak formulation is the fact that a conforming dis¬ 
cretization of the test space should be made by C^-elements. Hence Berggren proposed 
in [4] to introduce two new variables, (f := —Av and C ■= dnV allowing to perform a 
simpler numerical analysis, see Subsection 2.4. In Subsection 2.5 we propose another 
method that consists of regularizing the Dirichlet datum u by approximating it by a 
sequence of functions in 77^/^(F) using for example an interpolation operator. This 
allows to compute a sequence of weak solutions y^ G 77^ (D), and we show that they 
converge to the very weak solution with an explicit convergence rate. 

A negative result about the well-posedness of the weak formulation with L^(F)-data 
completes the discussion on the continuous level. 

Section 3 is devoted to the numerical analysis. We start with Berggren’s numeri¬ 
cal approach and recall his error estimates in Subsection 3.1 for completeness. Next 
we perform in Subsection 3.2 a numerical analysis of our regularization approach and 
prove error estimates for the piecewise linear approximation on a family of conforming, 
quasi-uniform finite element meshes. Notice that it turns out that on the discrete level 
Berggren’s approach is a particular case of our regularization strategy. The convergence 
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order is ^ in the convex case but smaller in the non-convex case. This reduction can be 
explained by the singular behaviour of the solution of the dual problem. In our paper 
[2] we investigate the singular complement method to remedy the suboptimality of the 
standard finite element method in non-convex domains. 

Finally, in Section 4 we present numerical tests in order to illustrate that our error 
estimates are sharp. The paper ends with some remarks about the three-dimensional 
case and about data with different regularity than assumed above. 


2 Analysis of the boundary value problem 

In this section we analyze several ways how to understand the solution of the boundary 
value problem (l.I) for which we cannot expect a weak solution y G 

For keeping the notation succinct we assume that the polygonal domain H has at most 
one non-convex corner with interior angle, denoted by oj. Let r, 6 be the corresponding 
polar coordinates and define \:= tt/oj. The boundary segments of are denoted by Fj, 
j = l,...,iV. 


2.1 Method of transposition in convex polygonal domains 

The method of transposition goes back at least to Lions and Magenes, [16], and is used 
by several other authors including French and King, [12], Casas and Raymond, [6], Dec- 
kelnick, Giinther, and Hinze, [11], and May, Rannacher, and Vexler, [17]. Since by partial 
integration the derivation 

= -{^y,v)n = (Vy, Vu)q for v E 

= {y,dnv)r - {y, Av)n for v E n H^{Q) 

is valid we get the very weak formulation-. Find 

yeL‘^{Q): {y, Av)n = {u,dnv)r - {f,v)n Vu E 77^(12) n iLg (12). (2.1) 

May, Rannacher, and Vexler, [17], proved the existence of a solution in L‘^{Q) in the case 
of a convex polygonal domain 12. 

Lemma 2.1. If the domain 12 is convex there exists a unique solution y E of 

problem (2.1) that satisfies the a priori error estimate 

provided that u E OjLi j)' ^ (^^(^) 72 ^7g(12))^ 

1 /2 

Recall that Hqq (Fj) is the space of functions whose extension by zero to F is in 
i7V2(r). 
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Proof. For being self-contained we sketch here the proof of [17, Lemma 2.1], The idea is 
first to assume more regular data, u G / G such that a weak solution 

y G exists, then to show (2.2), and finally to use a standard density argument since 

H\n) n 4 and FV2(r) 4 l^t) = nf=i 4 nf=i 

The estimate (2.2) is proven by using a duality argument. Due to the convexity of the 
domain there exists a solution w G 77^(0) of the auxiliary problem 

— Aw = y in D, tc = 0 on F, (2-3) 

such that 


\\y\\h(n) = -^w)n = if, w)n - {u, dnw)r 

nf.i4o"(ri)' + ll•^ll(^^2(o)n^^l(^))') 14142(0), 

where one uses that the mapping n Hq{Q) —> n^i^oo^(4)> ^ dnW, is 

surjective due to [13, Thm. 1.5.2.8]. The desired estimated (2.2) is then obtained by 
using the a priori estimate < c||7/11^2and by division by ||y|| 2 , 2 (Q). □ 

The method of proof of this lemma will be revisited in Section 2.5, where we start the 
discussion of the regularization approach. However, there we will work with different 
function spaces such that non-convex domains are included in the theory as well. The 
given proof of Lemma 2.1 is even restricted to convex domains since the isomorphism 

Au; G l 2(D), u;|r = 0 4^ w £ H^{n) n H^{n) 

is used. If the domain is non-convex one loses at least the uniqueness of the solution 
y of (2.1). For example, take H = {(rcos0,rsin0) G : 0 < r < 1,0 < 9 < lo} and 
A = tt/u!, then both yi = sin(A6l) and y 2 = r~^ sin(A0) are harmonic in D and yi = y 2 
on F. Both satisfy (2.1) with / = 0 and g = yi = y 2 on F. Hence, one needs a larger 
test space in order to rule out y 2 . 

2.2 Method of transposition in general polygonal domains 

In a hrst instance we replace the test space iL^(H) n F7o(H) by 

V := n Hi(n) with ffi(fi) -.= {u g H^(n) Av g (2.4) 

Since |u|/^i(q) < c\\Av\\l 2 (^^^, the graph norm in H, that is ||A7;||i2(Q) -b |u|hi(o), is 
equivalent to ||A7;||2,2 (q) such that we will use henceforth 

||u||y = ||A7;||i2(Q). 

Furthermore, let us denote by Vr the space of normal derivatives dnV of functions v G 
V. According to [13, Theorem 1.5.3.10] this space is well-defined and a subspace of 
F7QQ^(Fj)'. The natural norm in Hr is given by 

||c/||vy := inf{||7;||y :v £V,dnV = g} . 
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In particular, we have for v £ V that 


ll^nullvp < llulli/. (2.5) 

Since the previous definitions of the spaces V and Vp are rather formal let us discuss the 
structure of these spaces. 

Remark 2.2. The spaces V and Vr can be characterized as follows: 

1. If n is convex the spaces V and n Hq{Q) coincide. However, in non-convex 

domains the situation is different. In this case there is the splitting 

V = {H\n) n H^{n)) © Span{|(r)r^sin(A0)}, 

where ^ denotes a smooth cut-off function which is equal to one in the neighborhood 
of the non-convex corner. For more details we refer to [14, Sections 1.5, 2.3 and 
2.4] and [13, Theorem 4.4.3.7]. 

2. The mapping -£■ ^ dnW, is surjective due to [13, 

Thin. 1.5.2. 8]. Accordingly, in the convex case Hr is just whereas 

in the non-convex case there holds 

^r= ©Span{e(r)r"-i}. 

3. In the non-convex case there is H H^{Q)riHQ{Q) for s < 1 + A and Vr ^ 

for t < A — This implies n 77q(0))^ ^ V' and R*(r)' ^ 14p. 

Lemma 2.3. Let f £ V and ri G Vp. Then there exists a unique solution 

y £ L^(H) : (y, Av)n = {u, dnv)r - (/, v)n \/v £V (2.6) 

with 

\\y\\L‘2{n) < \\u\\v^ + WfWv- 


Proof. The proof of this lemma is based on the Babuska-Lax-Milgram theorem. Due 
to (2.5) the right hand side of (2.6) defines a linear functional on V. Moreover, the 
bilinear form is bounded on L^(D) x V. The inf-sup conditions are proved by using the 
isomorphism 

Av £ L^(D), Ujr = 0 v £ V. 

In particular, we obtain by taking y = Av 


|(j/,Au)n| ^ {Av,Av)n 
y&L‘2{n) \\y\\L'^{n) l|Au||r2('Q) 


l|Au||r2(t^) = Ijullv', 
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and by taking the solution v €V of Av = y 


Ky, An)n| ^ {y,y)n 
vev \\v\\v ~ \\y\\L^(n) 


\\y\\L2{n)- 


(2.7) 


The existence of the unique solution y G L^(n) of problem (2.6) follows from the standard 
Babuska-Lax-Milgram theorem, see for example [3, Theorem 2.1]. The a priori estimate 
follows from (2.7), and (2.6) and (2.5), 


hh^n) < sup 

v£V 


l(i/, A-;;)n| 
||u||y 


= sup 
v&V 


\{u,dnv)r - {f,v)n\ 
ll^llu 


< 


1^7 + 


V' 


□ 


Note that if a weak solution y G {u G : n|r = u} exists, then with the help of 

the Green formula 

{dnv,x)r = (Vu, Vx)n + {Av,xh Vu G Hi{n) D B, Vy G H\Q), (2.8) 

see Lemma 3.4 in the paper [9] by Costabel, it is also a very weak solution. (Set x = y 
and use (Vt>, Vy)n = {f,v)n.) 

A possible difficulty with this formulation is that a conforming discretization with a 
finite-dimensional space C V would require the use of G^-functions. This is simple for 
one-dimensional domains but requires a lot of degrees of freedom in two (and more) 
dimensions. 

We finish this subsection with a regularity result. 

Lemma 2.4. The unique solution y G of problem (2.6) with u G L^(r) and / = 0 

belongs to and to 

:= {z G L\n) : G (2.9) 

where d{x) is the distanee of x to the boundary T. Furthermore, there exists a positive 
eonstant c sueh that 

Il?/llj^i/2(r2) + ||'^^'^^Vy||i2(f^)2 < c||u||j;^2(r). (2-10) 

Remark 2.5. The book by Chabrowski, [7], deals exclusively with the Dirichlet problem 
with L? boundary data for elliptic linear equations. The solution is searched there in 
the Sobolev space (2.9) but domains of class were considered only. 

Proof of Lemma 2.f. In a first step, we use an integral representation and some prop¬ 
erties of the layer potentials to get a solution with the appropriate regularity. Theorem 
4.2 of Verchota’s paper [24] shows that the operator -\- K, K being the boundary 
double layer potential, is an isomorphism from L^(r) into itself (see also Corollary 4.5 
of [19]). According to the trace property for the double layer potential K, (see Section 1 
and Corollary 3.2 of [24]), there exists a unique harmonic function z such that 

z ^ u a.e. in nontangential cones 
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with the notation from [24], and is given by 

z = lC{y + K)-^u. (2.11) 

But due to the above mentioned isomorphism property and Theorem 1 of Costabel’s 
paper [9] (and its following Remark), we obtain that 

ll^llfl'i/2(n) ^ c\\u\\i2(^py (2-12) 

Note that Theorems 5.3, 5.4 and Corollary 5.5 of the paper [15] by Jerison and Kenig 
also yield 

11(5^/^V2;11^2(q) 2 < cllttllx, 2 (r). (2.13) 

The second step is to show that z is the very weak solution, hence by uniqueness, we 
will get y = z. The regularity of y and the estimate (2.10) follow from our first step. 
For that last purpose, we use a density argument. Indeed let Un G H^(^) be a sequence 
of functions such that 

Un ^ u in T^(r), as n —)> oo. (2-14) 

Consider Zn = + K)~^Un and let G L^(n) be the unique solution of (2.6) with 

boundary datum Un and right hand side / = 0 (that is in Then by the estimate 

(2.12) and Lemma 2.3, we get 

yn —t y in L^(n), as n —)■ oo, 

Zn ^ z in , as n —>• oo. 

Furthermore by Theorem 5.15 of [15] Zn satisfies 


■yZn = Un on F, 

where 7 is the trace operator from into Hence we directly deduce that 

Un = Zn and by the above convergence property we conclude that y = z. □ 

Corollary 2.6. The unique solution y G L^(H) of problem (2.6) with u G L^(r) and 
f G belongs to and to IT^’^(H) from (2.9). There exists a positive 

constant c such that 

\\y\\m/^in) + < c (llnllj;^2(r) + ||/||h-i(q)) • 

2.3 Method of transposition employing weighted Sobolev spaces 

Alternatively to the space n Hq{Q), one can use the test space 

F|’2(H)nHi(H), /3>1--, 

^ UJ 
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in the non-convex case, where is a weighted Sobolev space of the class 

:= j., e V'(n) : 


(2.15) 




where we use standard multi-index notation. For later use we also introduce 


L}{n) :=v^’\n). 

First derivatives of V^’^(0)-functions belong to by definition. The trace space 

of is n^i see [18, Lemma 1.2] or [20, Theorem 1.31]. In the next 

lemma, we will use the spaces 

for to > TT, 

L^(n) for U! < TT, 

2 2 9 

for /3 G (1 — 1]. We endow with the Vg’ (n)-norm for ui > tt and the Fr^(n)-norm 

for a; < vr, as well as Yp with the L^^(n)-norm for oj > tt and the L^(n)-norm for uj < tt. 

Remark 2.7. Let us discuss the dehnition of the space Vg and the restriction of the 
weight /3 to the interval (1 — 1] in the non-convex case: 

1. We require /3 G (1 — 1] in order to have the isomorphism 

Av G L^(n), U|r = 0 V € Vp. (2-16) 




•= J 


n Hq{Q) for u) > TT, 
n Hq{Q) for io < TT, 


2. It is possible to use the test space n Hq{Q) in convex domains as well. 

However, this implies a loss of information about the solution since this test space 
is smaller than H^{Q) n Hq{Q) due to the fact that the weight /? can be negative. 

3. Note that L^(Fj)for/3 < This implies L^(F) = ^]]][^;^ L^(Fj)^ 

• This means that L^-boundary data are included in the fol¬ 
lowing discussion if /? < 

Lemma 2.8. Let /? G (1 — 1] and assume that f £ Vp and u G (^O^i ■ 

Then there exists a unique solution 

y eYp^ L^{n) : {y, Au)o = {u, dnv)r - (/, v)n Vu G Vp (2.17) 


with 




< c 


V' 


- + Ill'll 
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Proof. The convex case was already treated in Lemma 2.1, hence we focus on the non- 
convex case. We proceed as in Lemma 2.3. The right hand side of (2.17) defines a 
continuous functional on Vg. Furthermore, the bilinear form is bounded on Yg x Vg, 

|(y,Au)o| < ||r"^y||L2(Q) ||r^Au||i2(Q) < ||y||y3 \\v\\vf,. 

The inf-sup conditions for the bilinear form are proved by using the isomorphism (2.16), 
in particular ||u||yg < c||Au|| 2 , 2 we obtain by taking y = r^^Au 

|(y,Au)Q| (r^Au,r^Au)n _ 

S \\y\\Yp “ iiAu|ii2(f2) - ^ ^ 

and by taking the solution u G Vg of Au = r~^^y with ||u||< c|| Au||£ 2 (f^^ = c||y ||^2 


sup 

veV/s 


\{y,^v)n 

Il'f^llV/s 


(r i^y,r l^y)n 

cWyhf^in) 




(2.18) 


The existence of the unique solution y G Lfi^iPL) of problem (2.17) follows now from 
the standard Babuska-Lax-Milgram theorem, see for example [3, Theorem 2.1]. The 
a priori estimate is obtained with (2.18) and (2.17), 


|y« < 


c sup 

v&Vp 


l(y, Au)o| 


’Wvp 


= c sup 

veV/j 


\{u,dnv)r - {f,v)n\ 


’Wvp 


This ends the proof since we already noticed that the enumerator defines a continuous 
functional on V)^. □ 


2.4 Berggren’s approach 

Berggren’s approach [4] avoids test functions in H^{Q) n Hq{0.) in an explicit way. It 
can be explained as if we substitute ip := —Av and C := dnV in (2.6), 

yeL'^in): {y,ip)Q =-{u,()r + {f,v)n 'iipeL^{Vl). (2.19) 

The relationship between ip G L?‘{Q) and both v £ V and C £ li" can be expressed by 
the weak formulation of the Poisson equation, 

veV: {Vv,V^I;)n = {p,f^h VV-G (2-20) 

and a reformulation of the Green formula (2.8) in the form 

CGPr: {C,x)r = {Vv,Vxh-ip,xh Vy G 77^(12) \ ^^(L!). ( 2 . 21 ) 

Note that Berggren’s formulation is not a system with three unknown functions since 
the second and third equations compute actions on the test function ip. Indeed, let 
S : L?{yi) —)• V and F : L^(Q) —>• Vr be the solution operators of (2.20) and (2.21), 
respectively, defined by Sip := v and Fip := (, then we could also write 

y G L^(I2) : (y, ip)n = -{u, Fip)r - (/, Sip)n \/ip G L^{n) 

instead of (2.19). 
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Lemma 2.9. Berggren’s formulation (2.19), (2.20), (2.21) is equivalent to the formula¬ 
tion (2.6). 

Proof. We first assume that y G LfiVt) satisfies (2.6) and show (2.19)-(2.21). For any 
G L^(rf) let V be the variational solution of — Au = g) defined by (2.20), hence v £ V 
and C := dnV £ Vr- Based on the formula (2.8), we obtain (2.21). With (2.6) we hnally 
get also (2.19). 

Let now y satisfy (2.19)-(2.21). Since G L?‘{Pl) we get from (2.20) that v £ V. 
Moreover, we obtain Q = dnV £ hr from (2.21). Hence equation (2.19) becomes 

-(y, Au)o =-(M,5nu)r + (/,u)o yv£V 

due to the isometry between L‘^{Q) and V. □ 

Remark 2.10. Berggren used the regularity v £ with some e > 0 which 

implies C, £ H^iT). However, he did not consider the maximal domain of the elliptic 
operator, i.e., v £ V and Q £ Vy- But with the explanations of Subsection 2.2 these 
regularities should be obvious. Thus the result of Lemma 2.9 is slightly more general 
than that of Berggren. 

2.5 The regularization approach 

A further idea is to regularize the boundary data and then to apply standard methods. 
This approach has already been considered within the proof of Lemma 2.1. In contrast, 
we do not use the isomorphism 

/\w£L^{n), w\y = Q ^ w £ 

which can only be employed in case of convex domains, but the isomorphism 

Au G L^(0), U|r = 0 V £V. 

This allows us to apply the regularization approach in the non-convex case as well. 
Moreover, we propose two different strategies how the regularized Dirichlet boundary 
data can be constructed in an explicit way. Thereby we will be able in Subsection 3.2 
to calculate approximate solutions of the regularized problems based on a finite element 
method. For the data we assume henceforth u £ L^(F) and / G This is not 

only for simplicity but also due to the fact that already for Dirichlet boundary data in 
L^(r) the convergence rates of the approximate solutions in Subsection 3.2 tend to zero 
as the maximal interior angle tends to 27r. 

We start with general convergence results for the regularized solutions. To this end 
let £ RV 2 (r) be a sequence of functions such that 

lim \\u — M^ll r 2 /r'i = 0 . 
h^o ^ ^ 

Let now y^ £ Yff := {u G H^{Pl) : u|r = u^} be the variational solution, 

y^£Y]f-. {Vy\Vv)n = {f,vh ( 2 . 22 ) 


10 


Lemma 2.11. Let u G L^(r) and f € H Then the limit y := lim exists, 

h—>-0 

belongs to LP‘{VL), and is the very weak solution, that means it satisfies (2.6). 

Proof. First we show that y^ is a Cauchy sequence in From (2.22) and Green’s 

formula, we have for any v € V, 


if, v)n = (Vy'*, Vv)n = -{y^, Av)n + {y^, dnv)r, 
if,vh = (Vy'*', Vu)o = -{y’^',Av)n + {y^',dnv)r. 

Hence due to y^ = and y^' = u^' on F, we deduce that 

[y^ - y ^', = [u^ - u^', dnv)r Vu G V. 

Now for any z G L^(H), let Uz G H be such that 

Avz = z, 


that clearly satisfies 

WdnVzhfir} < chzWH^in) < c ||2;|| l 2( s ^) 
with some s G (|, 1 + A), s < 2. Finally, we obtain with (2.23) and (2.25) 


(2.23) 


(2.24) 


(2.25) 


h'u {y -y ,z)n {u^ -u^ ,dnVz)r 

-y llL2(m = sup = sup — 

2eL2(o),2^o ll^llL2(n) 2eL2(n),2^o 


'IIl 2 (q) 


^11 h h'w ll^n'l’2||L2('p) ^ 

< ||u -u ||i 2 (r) sup ——2 -= c||u -u ||p 2 (r). 

2eL2(n),2^o IpIIl2(o) 


Since converges in L^(F), it is a Cauchy sequence and hence also y^ is a Cauchy 
sequence and converges in by the completeness of Lp‘{Ll). 

From V C 1/q(H) we obtain by (2.22) and the Green formula (2.8) 

(f,v}Q = (Vy^,Vv)Q = {Av,y^)n - {dnV,y^)r 
= (An, y'")n - {dnV, u^)t Vu G V. 

Since Av G LP‘{VL) and dnV G L^(F) we can pass to the limit and obtain that the limit 
function y satisfies (2.6). □ 

We can estimate the regularization error by a similar technique. 

Lemma 2.12. Let s = ^ifLlis convex and s € [0, A — ^) if Q is non-convex. Then the 
estimate 


\\y - y^llL 2 (n) < c\\u - 


holds. 
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Proof. We use the approach of the proof of Lemma 2.11 and just replace u^' by u and 
y^' by y to get 

{y-y^,l\v)n = {u-u^,dnv)T Vu G F. (2.26) 

Again, for any z G we let G 1^ be such that = z but estimate now in a 

sharper way 


||< 9 nU 2 ||//^( r ) < c\\Vz\\hs+3/2(^^-) < c\\z\\l2(^q) 
As in the previous proof we get 


(2.27) 


-y^llL 2 (n)= sup 


{u - U^,dnVz)i 


^ II II 

z£L'^{n),zj^0 ll^llL2(r2) 

, II fell l| 5 n^^ 2 :||H'*(r) 

< ||u - u llH-^(r) sup ——n- 

Zel'din),z^o \\z\\L‘2{n) 


< c\\u — u‘ 


h I 




Actually, the proof is for s > 0 but of course the statement holds when s is decreased. □ 


A choice for the construction of the regularized function could be the use of the 
L^(r)-projection lihU into a piecewise polynomial space on the boundary (which we call 
y)f in Section 3) or the use of the Carstensen interpolant ChU, see [5]. Namely, if Ar is 
the set of nodes of the triangulation on the boundary, we set 

ChU = ^ Tr^{u)X^, 
xsA/r 


where is the standard hat function related to x and 

fpuAx _ (u, Aa,)r 


7Tx(u) = 


fp Xx ( 1 , Xx ) r 


The advantages of the interpolant in comparison with the L^-projection are its local 
definition and the property 


u G [a, 5] => ChU G [o, 6], 

see [ 10 ]; a disadvantage for our application in optimal control is that ChUh 7 ^ Uh for 
piecewise linear Uh- We prove now regularization error estimates for the case that the 
regularized function is constructed via ChU or HhU. 

Lemma 2.13. If is the piecewise linear Carstensen interpolant of u or the L‘^{T)- 
projection of u into a space of piecewise linear functions, then there holds 

11^ “ '*^^lli7-^(r) < c/i^ll^llL2(r)) s G [0,1] 

as well as 

||tt^||L2(r) < c||tt||2,2(r). 
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Proof. The interpolation error estimate and stability result are derived in [5, 10] for 
domains. The proofs can be transferred to estimates on the boundary T. For the sake 
of completeness we sketch these proofs in the Appendix, see Lemma A.2. 

In the case of the L^-projection, the second estimate holds with constant one. For the 
first estimate, we notice that by using the properties of the L^-projection, we have 


{u - Uhu, (/?)r = {u- Uhu, ip - = {u,p- Iihp)v < IkllL^ir) \W - ^hp\\L'^{T) 

< \\u\\l'^{t)\\p - Chp\\L'^(T) < c/i^||it||L2(r)||y?||j^s(r) 


where we used Lemma A.2 in the last step. We conclude 


\u-u ||_H'-s(r) — 


sup 


{u - IihU,p)l 




< c/i*||u||i2(r). 


which is the assertion. 

By setting s = 0 in the previous lemma we obtain 

.h\ 


□ 


\u 


- u ||L2(r) < c||tt||i2(r) 


for the different choices of the regularized function u^. This means that for u G L^(F) 
the difference u — is uniformly bounded in L^(r) independent of h. However, we 
require strong convergence in L^(F) for u G L^(r), i.e. 

lim \\u — u^ll r 2 /r'i = 0. 
ft-s-o ^ ’ 

This is subject of the next lemma. A comparable result for the Ritz-projection can be 
found in e.g. [8, Theorem 3.2.3]. 

Lemma 2.14. Let u G L^(F) and let be the piecewise linear Carstensen interpolant 
of u or the L'^fT)-projection of u into a space of piecewise linear functions. Then there 
holds 

lim \\u — u^ll r 2 /r'i = 0. 
h^o ^ ’ 

Proof. We show the validity of this lemma for the Carstensen interpolant. The conver¬ 
gence result for the L^(r)-projection can be proven analogously. 

Due to the compact embedding Lf^(F) L^(r) there exists a sequence of functions 
Un G H^{T) such that 

lim \\u - UnWi^iD = 0 


with 


— ^ni 


(2.28) 


where the constant Cn may depend on n. Thus, for every e > 0 there is a positive integer 
N such that 

Ik —'*^n||L2(r) < - (2.29) 
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for n > N. By inserting the function Un and its Carstensen interpolant ChUn as inter¬ 
mediate functions into the desired term we obtain 

\\u - C'/in||i2(r) < ||n - UnWi^^r) + \\un - ChUn\\L^(r) + \\Chiun - u)\\L2{r)- 
The results of Lemma 2.13 and the inequalities (2.28) and (2.29) imply 

l|w - ChU\\L2^r) <c(^^+ hcn^ . 

Since for every constant Cn there is a parameter hn such that hcn < | for all h < hn we 
arrive at 

\\u - Chu\\L 2 (Y^ < ce 

and the desired result follows. □ 

2.6 A negative result 

Since the boundary datum u G L^(r) .b^(r) is trace of a function w G TT^’^(n) an 

idea could be to search y = yo + w with yo|r = 0 and 

yoGl<’^(L!): {Vyo,Vv)n = {f,vh - {Vw,Vv)n =: {F,vh yvGW^’°°{Q). 

However, bilinear forms a{u, v) : ITq^’^(H) x (H) —>■ C are investigated in [22, 23] 
and it is shown that 

VFg IT^’P'(H)* 3zG1To’^(H): aiz,v) = {F,v)n Vu G ITo’^'(H) 
holds for p G (1, oo) only since Cp —)• oo for p —?■ 1 in the inf-sup condition 

Cp sup > ||Vz|Up(o) for 1 < p < oo. 

3 Discretization of the boundary value problem 

Let {Th)hyQ be a family of conforming, quasi-uniform hnite element meshes, and intro¬ 
duce the hnite element spaces 

Yh = {vh^H\n)-.Vh\T^Vi'iT eTh], Yoh = Yhf^Hl{^), Y^ = Yh\an. (3.1) 

3.1 Berggren’s approach 


Let u G L^(r) and f € H ^(H). We discretize the formulation 

(2.19), (2.20), (2.21) in 

a straightforward 

manner, 


Vh^Yh'. 

{Uh-, ^h)n = -{u, Ch)r + if,Vh)n 

Vpft G Yh, 

Vh £ Y^h : 

(Vu/i,VV'h)Q = {^h,'il^h)n 

e Yoh, 

a e Y^ : 

{Ch,Xh)r = {Y^Vh,Vxh)n - {^h,Xh)n 

Vxh G Yh \ Yoh. 
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Note that Ch F^h and Vh / S(ph with F and S from Subsection 2.4, hence it is only 
an approximate Galerkin formulation. 

Berggren showed in Theorem 5.2 of [4] that this formulation is equivalent with the 
standard finite element approximation with L^-projection of the boundary data, 

Uh^Yh- {Yyh,'Viph)n = if,^h)n G lofe, (3.2) 

{yh,^h)r = {u,iph)r \/iph€Y^. (3.3) 

Note that yh = Uh on T when u = Uh & Yj^. This will be of interest in the discretization 
of optimal control problems. 

Berggren proved also the following discretization error estimate, see Theorem 5.5 of [4]. 

Lemma 3.1. Let to be the maximal interior angle of the domain Ll, and denote by 
X = tt/uj the corresponding singularity exponent. Let s' G (0, be a real number with 
s' < A — and let s G [0, s') be a further real number. Then the error estimate 

\\y — yh\\L^{n) < c + h^ lkllL2(r) + = 0{h^) 

holds, this means that we have convergence order s, 

. ) h ~ ^ convex domains, 

s — ixiiri 'I 2 ^ ^ 2 J ^ — I 1 

I A — ^ — e for non-convex domains, 

e > 0 arbitrary. 

Note that s —)■ 0 for w —>■ 27r. 

We will show in the next section that in the convex case the convergence order is 
without e. 

3.2 Regularization approach 

We consider a regularization strategy such that G T^, see Subsection 2.5. Recall that 
the corresponding solution y^ G Yff := {u G H^{LL) : u|r = u^'\ is defined via (2.22). For 
a regularization using the L^(r)-projection or the Carstensen interpolant we have the 
regularization error estimate 

\\y - 2 /^IIl 2 (o) < c\\u - u^Wn-^ir) < c/i^ll^i||L2(r) (3.4) 

with s = ^ if n is convex and s G [0,A — ^) if fl is non-convex, see Lemmas 2.12 and 
2.13. 

The finite element solution y^ is now searched in := n Y^ and is defined in the 
classical way, 


yh^Y^h- {Yyh,Vvh)fi = {f,Vh)n G Tq/i- (3.5) 

Note that, if we construct u^ by the L^-projection, we recover the Berggren approach 
as a special case. 
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Lemma 3.2. The finite element error estimate 


\\y^ - yh\\L2{n) < ch^ + lkllL2(r)) 

holds for s = ^ in the convex case and for s < X — ^ in the non-convex ease. 

Before we prove this lemma we can immediately imply the following final error estimate 
for this approach. By the triangle inequality we have 

\\y - y/illL2(o) < \\y - y’^\\L^(n) + l|y^ - yh\\L^(n)- 

The first term is already treated in (3.4). The second term is treated in Lemma 3.2. 
Corollary 3.3. The diseretization error estimate 

\\y - 2//i||l2(q) < ch"^ + ||n||i2(r)^ 

holds for s = ^ in the convex case and s G [0, A — j) in the non-convex case. 

Note that the order can be improved if the boundary datum u is more regular, see 
Remark 5.3. 

Proof of Lemma 3.2. Let Bhvfi G be the discrete harmonic extension defined by 

{VBhv!^, Vvh)n = 0 Vu/j G Toh 

which satisfies 

\\VBhu'^\\L2^n) < cWu’^Wny^r) < ch~^^‘^\\u^\\L2(r) < ch“^/^||u||i2(r). (3.6) 

The hrst estimate can be cited from [17, Lemma 3.2], the second follows from an inverse 
inequality, the third from Lemma 2.13. 

Now we notice that 

= BhU^ + as well as yn = BhU^ + yoh, 
where y(^ G and yoh G Toh satisfy 

{^y^,Vv)n = {f,v)n-{V{Bhu’^),Vv)n e H^{n), (3.7) 

(Vyo/i, Vu/j)o = if,Vh)n - {V{BhU^),Vvh)n G Tq/i- (3-8) 

Hence yoh is the Galerkin approximation of y^ and therefore (since Iq/i T Hq{LI)) 

(V(2/o - 2/o/i), Vuh)o = 0 Vu/j G Tb/i- (3.9) 

By taking u = t/q in (3.7) (resp. Vh = yoh in (3-8)), we see that 

11^2/0llL2(r2) ^ Il/lli^-i(o)ll2/o+ ll^(-®/*'“^)llL2(o)||Vyo||2,2(Q), 

W'^yohWhin) - ll/llH-i(n)l|yo/ilUi(o) + ||V(R/iU^)||i:,2(Q)||Vyoh||L2(Q). 
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By the Poincare inequality we obtain 


ll^ 2/0 IIl 2 (o) < c (^ 11 / 11 ^- 1 ( 0 ) + I|V(^/i^^)IIl 2 (o)^ , 
l|Vyo/i||L2(r2) < c (^||/||//-i(o) + ||V(i?/jn^)||i2(o)^ . 


With the help of (3.6) we arrive at 

ll^?/0 IIl 2(0) + l|Vy0fe||l,2(o) < c 

Now as before we start with 




+ h-V 2 | 




(i/o - i/O/i, z)q. 


lli/o - 2 / 0 fc||L 2 (O) = sup 

zeL2(o),2^o IpIIl2(o) 

Letting again G P be such that = z, we get 

(V(yo - yofe), V'^z)n 


lll/o-yoh||L 2 (o) = sup 

2 eL 2 (O), 27^0 

and therefore thanks to (3.9) we arrive at 


'I|l 2 (o) 


(3.10) 


-VOh\\L^(Q)= sup 

2 eL 2 (o), 2 ^o 


(V(i/o - Voh),'^{vz - hvz)h 

INIIl 2 (o) 


where Ih is the standard Lagrange interpolation operator. By the Cauchy-Schwarz in¬ 
equality, the well-known estimate 

\\vz - 4^^2||i,o < 'yz||//3/2+s(Q) < c/i^4+3||2:||j^2(o), 

with s = 2 in the convex case and s < A — ^ in the non-convex case, and the estimate 
(2.27) we obtain 

ho - 2 / 0 /i||l 2 (o) < ch^ 4 + 3 ||v(y^ - yoh)\\mn)- 
Using the a priori estimate (3.10), we conclude that 

ho ~ yoh\\L‘2{n) < (^^'^^11/11^-1(0) + ll^llL 2 (r)) • 


□ 

Note that we have now proved the convergence order ^ for Berggren’s approach, in 
the convex case. 
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Solution of the Problem 



Figure 1: Visualization of the exact solution 


4 Numerical test 

This section is devoted to the numerical verihcation of our theoretical results. For 
that purpose we present an example with known solution. Furthermore, to examine 
the influence of the corner singularities, we consider different polygonal domains 
depending on an interior angle a; G (0, 27r). These domains are defined by 

:=(-!, 1)^ n {x G : (r(x), (p{x)) G (0, V2] x [0, w]}, 

where r and tp stand for the polar coordinates located at the origin. The boundary of 
is denoted by F^ which is decomposed into straight line segments Tj, j = 1,... 
counting counterclockwise beginning at the origin. As numerical example we consider 
the problem 

-Ay = 0 in 

y = u onTj, j = 1 ,..., m(a;). 

The boundary datum u is chosen as 

u := sin(—0.4999(/9) on F^^. 

This function belongs to L^(F) for every p < 2.0004. The exact solution of our problem 
is simply 

y = sin(-0.4999(/?), 

since y is harmonic. A plot of it can be seen in Figure 1. We solve the problem numer¬ 
ically by using a finite element method with piecewise linear finite elements combined 
with either the L^-projection or the Carstensen interpolant of the data on the boundary. 
The finite element meshes for the calculations are generated by using a newest vertex 
bisection algorithm as described in [21]. The discretization errors for different mesh sizes 
and the experimental orders of convergence are given in Tables 1-3 below for interior an¬ 
gles io G {37r/4, 37r/2, 3557r/180}, where the discrete solutions based on the L^-projection 


18 



Solution of the Problem 


Solution of the Problem 



Figure 2: Visualization of the approximate solutions with /i = 1/8 


and the Carstensen interpolant are denoted by yh ,2 or yh,Cj respectively. We note that 
the errors are calculated by an adaptive quadrature formula. Apparently, the results are 
very much in congruence with the predicted orders. The different numerical solutions 
yh ^2 and yh^c for w = 37r/2 and /i = 1/8 are displayed in Figure 2a and Figure 2b, 
respectively. We see that the infinite boundary value in the origin is replaced by a finite 
one. If using the L^-projection of the data, we see also that the zero boundary values at 
the edge with 0 = 0 are replaced by an oscillating function which is typical for this kind 
of regularization. By using the Carstensen interpolant as regularization of the data, this 
can be avoided according to the local definition of this interpolation operator. But we 
note that these oscillations are a feature of the former regularization approach and do 
not disturb the approximation order. 


5 Extensions 

Remark 5.1. The paper is written for two-dimensional domains. However most results 
also hold in the three-dimensional case or can simply be extended to this one. The 
most crucial issue is the regularity which we have used for the corresponding adjoint 
problem since in three dimensional domains not only corner singularities but also edge 
singularities need to be taken into account. 

Remark 5.2. We assume / G for the discretization error estimates of Corollary 

3.3. This is only for simplicity. For defining the very weak solution and the numerical 
solution we only need / G {Hq{Q) n H^{Q)y n such that we could admit right hand 
sides from L^{Q) and Dirac measures as well. However, in this case the discretization 
error analysis demands an adapted proof, which exceeds the scope of this paper. 

Remark 5.3. We assume u G T^(F) for simplicity. The case u G R*(F), t > 0, is 
also of interest in the analysis of Dirichlet control problems, see [1]. The results can be 
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mesh size h 

^ unknowns 

\\y 

- L2(q^) 

eoc 

\\y 

- y/i,c l2(q^) 

eoc 

0.50000 

19 


0.26142 



0.26794 


0.25000 

61 


0.18577 

0.49289 


0.18973 

0.49794 

0.12500 

217 


0.13172 

0.49600 


0.13426 

0.49899 

0.06250 

817 


0.09331 

0.49745 


0.09497 

0.49940 

0.03125 

3169 


0.06605 

0.49838 


0.06717 

0.49965 

0.01562 

12481 


0.04674 

0.49902 


0.04750 

0.49982 

0.00781 

49537 


0.03306 

0.49942 


0.03359 

0.49992 

0.00390 

197377 


0.02338 

0.49967 


0.02375 

0.49998 

expected 




0.5 



0.5 

Table 1: Discretization errors and experimental orders of convergence (eoc) for 

w = 37r/4. 

Expected convergence rate: 1/2 





mesh size h 

^ unknowns 

\\y 


eoc 

l|y 

- yh,c\\L^(n^) 

eoc 

0.50000 

33 


0.73622 



0.77007 


0.25000 

113 


0.64484 

0.19118 


0.67086 

0.19897 

0.12500 

417 


0.56841 

0.18201 


0.58915 

0.18737 

0.06250 

1601 


0.50328 

0.17555 


0.52022 

0.17950 

0.03125 

6273 


0.44674 

0.17194 


0.46091 

0.17464 

0.01562 

24833 


0.39711 

0.16987 


0.40920 

0.17166 

0.00781 

98817 


0.35330 

0.16865 


0.36376 

0.16982 

0.00390 

394241 


0.31448 

0.16793 


0.32362 

0.16868 

expected 




0.16667 



0.16667 

Table 2: Discretization errors and experimental orders of convergence (eoc) for 

u) = 37r/2. 

Expected convergence rate: 1/6 





mesh size h 

^ unknowns 

\\y 

“ 2//i, 2 L2(f7^) 

eoc 

\\y 

- 2/fe,c l2(q^) 

eoc 

0.50000 

46 


1.1049 



1.1141 


0.25000 

159 


1.0693 

0.04721 


1.0732 

0.05406 

0.12500 

589 


1.0491 

0.02749 


1.0513 

0.02967 

0.06250 

2265 


1.0367 

0.01715 


1.0384 

0.01782 

0.03125 

8881 


1.0281 

0.01207 


1.0296 

0.01226 

0.01562 

35169 


1.0213 

0.00956 


1.0228 

0.00962 

0.00781 

139969 


1.0154 

0.00832 


1.0169 

0.00834 

0.00390 

558465 


1.0100 

0.00771 


1.0114 

0.00772 

expected 




0.00704 



0.00704 


Table 3: Discretization errors and experimental orders of convergence (eoc) for u = 
3557r/180, Expected convergence rate: 180/355 — 1/2 = 1/142 ps 0.007 
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improved in this more regular case: With s from Corollary 3.3 we have 

\\y - y^h^in) < ch - < ch^^^\\u\\Ht(r), t G [0, i]. 

The first step is again the application of Lemma 2.12, while the second can be proved 
in analogy to Lemma 2.13. The necessary prerequisites are already provided in Lemma 
A.2 and Remark A.3. Furthermore, when we check the proof of Lemma 3.2 we find that 
we obtain 

Wvo - yoh\\L^(n) < c/i® ^ £ [0, 5]- 

Hence it remains to prove 

||M^||//t(r) < c\\u\\Ht(r), (5-1) 

in order to conclude 

\\y — yh\\L‘ 2 {n) < *ll/lli^-i(f2) + Iklli^qr)) ) 

The estimate (5.1) is known for domains and can be proved for t G [0,1] also for the 
boundary T by using the inverse inequality and the approximation properties of and 
the Ritz projection. To this end recall that Tj, j = 1,..., A^, are the boundary segments 
of T, and let be the Ritz projection on Tjf Ir^ . Then we have 


\u 


N N 

^ 11 “^ “ WPi'^Wn^iTj) 

i=i i=i 

N N 

< ch-^ \\u'^ - Plu\\L2(r,) + Y ll'“II^^MR) 

i=i i=i 

N 

< ch~^\\u - u^h^ir) + Y 11^ “ Ph'^Wpirj) + IkllHqr) 

i=i 




i. e., we have (5.1) for t = 1. Since we proved (5.1) for t = 0 in Lemma 2.13 we get the 
desired result by interpolation in Sobolev spaces. 


Appendix: Error estimates for the Carstensen interpolant 


Recall from Subsection 2.5 that the piecewise linear Carstensen interpolant is defined 
via 


ChU 


Y T^x{u)\x-, T^x{u) 
x^Jsfr 


Jp 

fr 


(u, A3;)r 
(I5 ^a;)r 


where Ax is the standard hat function related to x. 
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Lemma A.l. The piecewise linear Carstensen interpolant satisfies the error estimate 

\\ip - ChTWi'^ir) < ch^\\T\\H’>{r), (A.2) 

foripeHfiV), se[0,i]. 

Proof. The interpolation error estimate (A.2) is in principle contained in [5], however, 
it is there an estimate on a domain such that we sketch the proof here. 

For s = 0 this estimate follows from the stability property 

hxiu)\\L2(uj^) < c||«||l2(^^), (A.3) 

where iOx is the support of on F. 

For s = 1 we use that 'Kx{w) = w for all constants w such that 

||w - vra:(t6)||i2(^^) = ||(n - u;) - -Kxiu - u;)||l 2(^^) < c\\u - U^||l 2(^^) < ch\\u\\Hi{y:^) 

via the Deny-Lions lemma, see also [10, Lemma 4.3] where the piecewise affine and 
Lipschitz continuous transformation of Ux to some reference domain is discussed in detail. 
For a boundary edge with end points xi and X 2 we have 

2 

IIW - C'/x^^lll,2(e) = ll(« - '^xAu))Xxi + {U- TTx2{u))>^X2\\L2{e) < ||w - VTa;, (n) || i2(e) . 

i=\ 

From these two estimates we obtain (A.2) in the case s = 1. In the remaining case 
s G (0,1) the error estimate (A.2) follows by interpolation of Sobolev spaces. □ 

Lemma A.2. The piecewise linear Carstensen interpolant satisfies the error estimate 

\W-Chu\\H-s{r) <ch^^^\\u\\Ht(T), 'SG[0,1], fG[0,1], 

as well as the stability estimate 


\\Chu\\L2(r) < c||tt||x,2(r). 

Proof. The second estimate follows directly from the fact that 0 < < 1 and the 

stability property (A.3). 

The first estimate is in principle contained in [10], however, it is there an estimate on 
a domain such that we sketch the proof here. First we notice that the definition of 7rx{u) 
is equivalent to (u — TTxiu), Aa;)r = 0 and hence we have 

{u-TTx{u),7rxiip)Xx)r = 0 G H%T). 

With this identity and with A^, = 1 we get 

(u - ChU,ip)r = {u'^Xx -'^7rxiu)Xx,^p)r = '^{u - 7rxiu),ipXx)r 

XX X 

= - TTx{u), {(p - Tix{p))Xx)r 

X 

< ^ ||u - vr^(n)||i2(^^)||<^ - ■nx{p)\\Lfiuj,) 

X 
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where ojx is again the support of on F. With similar arguments as in the proof of 
Lemma A.l we conclude 


{u - Chu,ip)r < c/i*+*||u||H‘(r)ll<FllH-(r)- 
By the definition of the negative norm, 

(u - ChU, ip)r 


\u-Chu\\H-s^Y') = sup 




we obtain the assertion of the lemma. 


Remark A.3. Note that this error estimate holds also for the L^-projection, 
||u - < c/i^+‘||u||Ht(r), sG[ 0, 1], tG[0,1]. 

It can be proved similarly by using (u — HhU, (p)r = {u — HhU, ip — n/i(/9)r. 


□ 
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